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Introduction 

Entry  of  hydrogen  into  metals  is  of  serious  concern  to  metallurgists  and  engineers,  since  it 
severely  degrades  the  mechanical  properties  of  metals  (1).  These  problems  arise  in  the  cathodic 
protection  of  metals,  power  plants,  and  in  environments  where  H2S  is  present  as  in  petroleum  refining 
(2)  where  the  hydrogen  evolution  reaction  (h.ex)  and  hydrogen  permeation  reaction  take  place  in  the 
corroding  or  cathodically  polarized  metal.  By  performing  hydrogen  charging  experiments  of  thin 
samples  using  the  Devanathan-Stachurski  cell,  the  permeation  characteristics  have  been  extensively 
studied  (3,4).  The  present  paper  seeks  to  analyze  the  h.e.r.  mechanism  and  to  predict  the  relationship 
between  the  permeation  flux  and  the  charging  and  evolution  (recombination)  fluxes.  A  thorough 
development  of  the  model  and  actual  computations  of  rate  constants  and  hydrogen  coverages  will 
appear  elsewhere  (5). 

Analysis 

Essentially,  three  steps  are  involved  during  cathodic  hydrogen  charging  of  metals.  They  are:  (1) 
the  hydrogen  discharge  reacdon  (proton  tunneling),  (2)  hydrogen  recombination  reaction  either  by 
chemical  recombination  or  electrochemical  desorption,  and  (3)  hydrogen  permeation  (mainly  by  bulk 
diffusion).  Of  these  three,  the  bulk  diffusion  step  is  usually  the  slowest .  A  detailed  schematic  of  the 
reactions  is  given  in  Fig.  1  below. 

The  charging  current  (ic)  is  given  by 

ic-io'd-Os)  e’aotT1  (1) 

The  hydrogen  evolution  current  (if)  (assuming  chemical  recombination  of  H  atoms)  is  given  by 

ir  =  Fk3  0s2  (2) 
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The  steady  state  hydrogen  permeation  current  (ioo)  is  given  by 

i«,  3  F  ~  Cj  (3) 

In  equations  1, 2  and  3:  io'  ■  Fkj  « io/(l-8e);  k>  ■  *e  exchange  current  density;  6c  -  the  equilibrium 
surface  coverage  of  hydrogen ;  ki»  the  discharge  rate  coefficient  ■  k°  c^+e  -aotE«l  (acid)  * 

k°e*®®Ecq  (alkaline);  k°  ■  the  rate  constant  for  the  forward  reaction;  cjj+  =  H+  ion  concentration;  a  = 

F/RT  -  38.94  (volts)*1  at  T  -  300K;  a  -  the  transfer  coefficient;  -  the  equilibrium  potential  for  the 

h.ex,  0s  » the  surface  coverage  of  hydrogen;  T|  -  the  hydrogen  overvoltage  »  E  applied  -  E°fc  ^3  =  the 
recombination  rate  coefficient;  Di « the  hydrogen  diffusion  coefficient  in  die  metal;  L  ■  the  membrane 
thickness. 

The  model  considers  a  selvedge  reaction  (as  a  result  of  proton  tunneling)  that  is  quite  fast  and 
constitutes  a  transition  layer  of  a  thickness  that  could  range  upward  from  ~  1  nm  (that  has  yet  to  be 
determined  by  special  experiments).  This  establishes  a  metal  subsurface  hydrogen  concentration,  cj,  at 
the  boundary  of  the  selvedge.  Thus,  hydrogen  diffuses  out  to  either  surface,  though  mostly  to  the 
charging  side  to  recombine  to  form  H2  molecules.  An  equilibrium  will  be  established  between  the 
surface  covered  (adsorbed)  hydrogen  atoms  and  hydrogen  just  below  the  surface  (in  the  adsorbed  state, 

with  concentration  Cs),  This  equilibrium  has  been  analyzed  before  (6,7)  giving  0s  =  Cs/k',  where  k’  = 
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the  equilibrium  absorption  -  adsorption  constant  and  cs  ■  cj  -  Cg.  (However,  it  is  to  be  emphasized  here 
that  the  selvedge  reaction  is  not  critical  to  the  development  and  application  of  this  model  Consideration 
of  selvedge,  on  the  other  hand,  helps  to  generalize  the  model;  in  the  absence  of  a  selvedge,  eg  ■  0  and 
c**ci).  Using  these  relations  along  with  equations  (1),  (2)  and  (3)  one  can  strive  at  the  following  set 
of  equations,  which  for  the  first  time  have  taken  into  account  die  effect  of  w  on  the  h.ex  kinetics: 


where  b  ■  L/(FDi)  ■  a  constant  for  a  metal;  also, 

1 120  mv  1  1-  dn  .1  !  *1  I  ,  I  jn 

120d^-  disrd <  diSFri <2  dioa 


!d^<!d^<^jd^<r240^ 


for  the  model. 


The  transfer  coefficient  is  given  by 


Idlnic) 

os-  hr/a 


i  i 


Details  of  the  derivations  are  shown  elsewhere  (3). 

The  potential  range  (qj.  to  r^)  in  which  the  recombination  reaction  will  be  coupled  with  the 
discharge  reaction  (commonly  observed  on  many  metals)  is  given  by  (5) 

T|c  -  [in  (iOkj/kj]  i  (ao)  (7a) 

«  Qn  (kjAlOkj))]  /  (aa)  (7b) 

where  the  superscripts  1  and  u  refer  to  the  lower  and  upper  limits  of  the  overpotential  range. 


In  hydrogen  permeation  experiments,  ic  is  set  and  when  the  permeation  current  becomes 
independent  of  time,  ioo  is  measured.  Then  if  =  ic  -  w 

If  the  plots  of  ioo  vs  Vir  (equation  (4))  and  ic  e  a(XTl  vs  [ioo  -  (equation  (5))  are  linear,  then  all 

of  the  coefficients  k',  k3t  Cg  and  io'  can  be  calculated.  Such  calculations  have  been  done  on 
experimental  data  from  the  literature  for  iron  and  nickel  membranes  and  found  to  verify  the  model  (5) 
since  these  two  plots  are  linear.  An  example  of  such  an  analysis  is  given  for  the  polarization  and 
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permeation  Haw  of  Bockris  et  al  (3)  obtained  on  Aimco  iron  in  0.1N  H2SO4  solution.  Fig.  2  shows 

the  plot  of  1m  vs  Vir  and  Fig.  3  shows  ic  «  v*  [i~  -  plotfor  die  data  of  Bockris  etal  (3).  It  is 
easily  seen  that  these  plots  are  linear  and  hence  the  model  can  be  applied  to  detennine  the  rate  constants 
and  exchange  current  density  (Table  1).  Then,  from  Eqn.  (2),  6s  can  be  calculated  using  the  k'  and  io 
values  obtained  from  the  slope  and  intercept,  respectively,  of  Fig.  3,  die  k3  value  obtained  from  the 

slope  of  Fig.  2  and  the  ir  value  from  ic  -  w  The  surface  coverage  (6$)  vs  the  hydrogen  overvoltage 

(T|)  plot  (for  the  Bockris  et  al  data)  is  shown  in  Fig.  4.  It  can  be  seen  that  the  coverage  is  quite  low  in 
the  potential  range  of  experimentation. 


7 


J 


FIG.  2.  Analysis  of  Data  of  Bockris  et  al  (3). 


FIG.  3.  Analysis  of  Data  of  Bockris  etal  (3). 


FIG.  4.  H  Coverage  for  the  Data  of  Bockris  et  al  (3). 


TABLE  1 

Calculated  Values  by  Applying  the 
Model  to  the  Data  of  Bockris  et  al  (3) 

io  *  0.5  (lA/sq.cm. 
ki  »5.5x  lO’^mol  /  cm2.s 
k3  ■  2.3  x  10"3  mol  /  (cm^.s 
k’»2.6x  10'5  mol  /cm3 

nj.*  -600 mV;  ti“«  -810mV 
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Thus,  this  model  effectively  accounts  for  the  contribution  of  ieo  to  the  overall  kinetics.  The  most 
important  prediction  of  this  model  is  that  io*  is  proportional  to  Vif  and  not  to  ■'Ac  as  assumed  in  earlier 
models  (3).  This  type  of  relationship  (io*  a  Vir)  has  been  previously  observed  (8). 

The  above  relationships  assume  that  q  »  RT/F  so  that  the  backward  reactions  can  be  neglected. 
Also,  the  Langmuir  isotherm  of  hydrogen  coverage  was  utilized  in  order  to  simplify  the  derivation. 
However,  in  many  cases,  the  reactions  (discharge  and  recombination)  are  activated  in  which  case 
Fntmlrin.Temlrin  corrections  (9)  have  to  be  applied  for  0j  in  the  equations  for  ic  and  ir.  Equation  (1) 
then  becomes: 

4  -  4(1-0.)  e^'e-0”  (8) 

and  equation  (2)  becomes 

4  ■  Fk30ge2a/e*  ^ 

where  /  ■  yRT,  y  being  the  gradient  of  the  apparent  standard  free  energy  of  adsorption  with  coverage. 
The  value  of  /  *  4  to  5  for  H  coverages  (9).  In  the  problem  of  enhanced  hydrogen  entry  in  the 
presence  of  H2S,  such  considerations  have  been  shown  to  be  necessary  (10).  The  modified 

relationships  between  io*  ir  and  ic  are  given  by 


and 

where 


In 


Jk) 


— J  i„  -  In 


k‘ 


In  (/(ic  ij)  *  o  a  (-T|)  +  In  (4') 


/(i*»U 


(afl>U 

k* 


a-f=> 


(10) 

(11) 


Thus,  equation  (10)  tells  us  that  i«o  will  not  be  lineariy  related  to  Vir  when  />o,  meaning  the  discharge 
and  recombination  reactions  are  activated,  probably  due  to  a  side  reaction  of  H2S  with  a  hydrated 
electron,  e‘aq  (10),  as  follows: 


and 


A 


I^S  +  ei,  -*  HjS* 
H2S'  +  H+  +  M  -»  M-H  +  H2S 


(12) 

(13) 


If  In  (- — )  vs  io»  and  In  (/(ici«.))  vs  q  are  linear,  then  this  side  reaction  and  overall  mechanism  can  be 

loo 


said  to  be  operating.  Then,  the  coefficients  k',  k3,  a  and  io'  can  be  computed  from  the  slopes  and 
intercepts  of  these  plots,  in  conjunction  with  the  iterative  solution  of  equations  (10)  and  (1 1). 
Eventually,  0S  vs  q  can  be  plotted;  and  0e,  io  and  kj  can  be  computed. 

Once  again,  the  potential  ranges,  where  the  recombination  and  discharge  reactions  are  coupled, 
can  be  estimated  from  equations  (7a)  and  (7b).  With  increasing  H2S  concentration,  these  potential 
ranges  have  been  found  from  analysis  of  available  data  according  to  the  above  model  to  become  less 
negative,  k3  progressively  decreases  suggesting  a  decreasing  surface  diffusivity  of  Had  atoms,  and  ki 
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progressively  increases  suggesting  that  the  discharge  reaction  is  enhanced  by  die  side  reaction,  equation 
(12),  followed  by  equation  (13).  This  will  also  explain  why  the  overvoltage  actually  decreases  rather 
than  increases  in  the  presence  of  H2S. 


Conclusions 


*  !n 


A  summary  of  a  recently  completed  analysis  of  hydrogen  electrode  reactions  during  aqueous 
cathodic  charging  of  metals  is  presented.  The  analysis  for  tne  first  time  takes  into  account  the  effect  of 
hydrogen  permeation  into  the  metal  on  the  h.ej.  From  the  model  all  of  the  kinetic  parameters  are 
computable  without  use  of  any  adjustable  parameters.  For  die  first  rime,  surface  coverages  and  rate 
constants  are  determinable  from  the  measured  charging  and  permeation  currents.  The  enhancement  of 
hydrogen  entry  in  the  presence  of  poisons,  such  as/HjSrcan  be  analyzed  with  the  model,  taking  into 
account  the  Fntmkin-Temkin  isotherms  in  the  discharge  and  recombination  reaction  kinetics.  ^  /  '  • 

Vtfei  rt  '*  r ,  •  i  V'r  ■*  •  J  t  i  ' 
(  Ju  '< !>  <  i  1  j  I  f y.. 
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